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GLOBAL EXISTENCE OF SOLUTIONS TO REACTION DIFFUSION 
SYSTEMS WITH MASS TRANSPORT TYPE BOUNDARY CONDITIONS 


VANDANA SHARMAt AND JEFF MORGAN* 

Abstract. We consider a reaction-diffusion system where some components react and diffuse on the bound¬ 
ary of a region, while other components diffuse in the interior and react with those on the boundary through 
mass transport. We establish local well-posedness and global existence of solutions for these systems using 
classical potential theory and linear estimates for initial boundary value problems. 
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1. Introduction. The idea that reaction-diffusion phenomena is essential to the growth 
of living organisms seems quite intuitive. Indeed, it would be rather hard to envision how any 
organism could grow and operate without moving its constituents around and using them in 
various bio-chemical reactions [16]. For example, bacterial cytokinesis is one process which 
can be modeled by reaction-diffusion systems. During the bacterial cytokinesis process, a 
proteinaceous contractile ring assembles in the middle of the cell. The ring tethers to the 
membrane and contracts to form daughter cells; that is, the “cell divides”. One mechanism 
that centers the ring involves the pole-to-pole oscillation of proteins Min C, Min D and Min E. 
Oscillations cause the average concentration of Min C, an inhibitor of the ring assembly, to be 
lowest at the midcell and highest near the poles [35], [27]. This centering mechanism, relating 
molecular-level interactions to supra-molecular ring positioning can be modelled as a system 
of semilinear parabolic equations. The multi-dimensional version of the evolution of the Min 
concentrations can be described as a special case of the reaction-diffusion system 
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where D is a bounded domain in R", n > 2, with smooth boundary M, A and Am denote the 
Laplace and Laplace Beltrami operators, 77 is the unit outward normal vector to D at points 
on M, and D and D are k x k and m x m diagonal matrices with positive diagonal entries 
{dj}i<j<k and respectively. F : x R"* —>• R"*, G : R^ x R"* —>• R^, iL : R^ —>• R^, 

and uq and vq are componentwise nonnegative smooth functions that satisfy the compatibility 
condition 


= G{uo, 'Po) on M. 
or] 
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For this model, n may represent the cell cytoplasm and M may represent its membrane. There 
are some components that are bound to the membrane, and other components that move freely 
in the cytoplasm. Also, the components on the membrane and cytoplasm react together on 
the membrane through mass action and boundary transport. In Section 7, we present two 
applications associated with (1.1), with one modeling the chemical reaction involving Min 
protiens for positioning of the ring, explained in [27]. We point out the study in [35] that also 
modeled these reactions. 

In general, system (1.1) is somewhat reminiscent of two component systems where both of 
the unknowns react and diffuse inside II, with various homogeneous boundary conditions and 
nonnegative initial data. In that setting, global well-possedness and uniform boundedness has 
been studied by many researchers, and we refer the interested reader to the excellent survey of 
Pierre [23]. 

In the remainder of the introduction, we assume H = 0 and k = m = 1. A fundamental 
mathematical question concerning global existence for (1.1) asks, what conditions on F and 
G will guarantee that (1.1) has global solutions, and how are these conditions related to the 
results listed in [23]? The focus of this paper is to give a partial answer to this question and to 
apply our results to (1.1). 

From a physical standpoint, it is natural to ask under what conditions the solutions of (1.1) 
are nonnegative, and the total mass is either conserved or reduced. It is also important to ask 
whether these conditions arise in problems similar to the above mentioned cell biology system. 
Conditions that are similar in spirit to those given in [20], [13] and [23] result in nonnegative 
solutions for system (1.1). More precisely, (1.1) has nonnegative solutions for all choices of 
nonnegative initial data uq and vq if and only if F, G, and FI are quasi-positive. That is 
F(a, 0), 0(0, a) > 0 whenever a > 0 (recall from above that 77 = 0 in the remainder of this 
introduction). Also, some control of total mass can be achieved by assuming there exists a > 0 
such that 

(1.2) F{u,v) + G{u,v) < a{u + V + 1) forallu,?;>0. 

Assumption (1.2) (discussed later), generalizes mass conservation by implying that total mass, 
u(x, t) dx + u(^, t) da, grows at most exponentially in time t. 

We suspect that the natural conditions, quasipositivity and conservation of mass, are not 
sufficient to obtain global existence in (1.1), and that it is possible to construct an example 
along the same lines as constructed in [24]. To this end, we impose a condition similar to 
Morgan’s intermediate sums [21] and [22]. Namely, there exists a constant Kg > 0 such that 

G(C, v) < + 1 ^ + 1) for all C > 0- 

In addition, we adopt a natural assumption of polynomial growth, which has been considered 
in the context of chemical and biological modeling (see Horn and Jackson [14]). That is, there 
exists I G N and Kf > 0 such that 

F{u,v) < Kf{u + V + ly for all v>0,u>0. 

In our analysis, we extend recent results of Huisken and Polden [25], Polden [15], and Sharpies 
[31] associated with x (0,T)) results for solutions to linear Cauchy problems on a 

membrane. We also verify and make use of a remark of Brown [4] which states that if d > 0 
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and the Neumann data 7 lies in Lp{M x (0, T)) for p > n + 1, then the solution to 



(fit = dA(p 

X £ id, 

0<t<T 

(1.3) 

dip 

X £ M, 

0<t<T 


p = 0 

X £ il, 

t = 0 


is Holder continuous on H x (0,T). We provide the proof of this result in section 5 for com¬ 
pleteness of our arguments. 

Note that the results of Amann [3] can be used to guarantee the local well posedness of 
(1.1) subject to appropriate conditions on initial data and on the functions F and G. However, 
those results do not provide the explicit estimates that are needed in our setting. Our approach 
keeps the analysis on comparatively simpler Lp spaces. 

It is worth mentioning that some of the results in section 5 are valid for domains that are 
only . Handling cases with weak smoothness conditions on curves or domain boundaries was 
one of the motivations for the results obtained in [4], [5], [9] and [10] , and these results may be 
of independent interest. 

2. Notations, Definitions and Preliminary Estimates. Throughout this paper, n > 
2 and H is a bounded domain in R" with smooth boundary M (9f2) belonging to the class 
with fi > 0 such that H lies locally on one side of its boundary, rj is the unit outward 
normal (from O) to M, and A and Am are the Laplace and the Laplace Beltrami operators, 
respectively. For more details, see Rosenberg [30] and Taylor [34]. In addition, m,k,n,i and j 
are positive integers, D and D are k x k and m x m diagonal matrices with positive diagonal 
entries {dj}i<j<k and {di}i<i<m, respectively. 

2.1. Basic Function Spaces. Let H be a bounded domain on with smooth boundary 
such that B lies locally on one side of dB. We define all function spaces on B and Bt — Bx{ 0,T). 
Lq{B) is the Banach space consisting of all measurable functions on B that are q^^{q > 1) power 
summable on B. The norm is defined as 

\\u\\q,B = \u{x)\‘>dx^ 

Also, 

Halloo,s = esssup{lit(x)l : x £ B}. 

Measurability and summability are to be understood everywhere in the sense of Lebesgue. 

If p > 1, then Wp(B) is the Sobolev space of functions u : H —>■ R with generalized deriva¬ 
tives, (in the sense of distributions) [s] < 2 belonging to Lp{B). Here s = {si, S 2 ,---,Sn), |s| = 
Si -b S 2 + .. + Sn, [s] < 2, and 9® = where di = The norm in this space is 

2 

ii«iiS= E ii^>iu 

ld=o 

Similarly, Wp’^(Bt) is the Sobolev space of functions u : Ht —R with generalized deriva¬ 
tives, d^d[u (in the sense of distributions) where 2 r -b [s] < 2 and each derivative belonging to 
Lp(Bt)- The norm in this space is 

2 

Mpk= E 

2r+|s|—0 
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In addition to we will encounter other spaces with two different ratios of upper 

indices, W 2 ^’°(St), W 2 ^{Bt), V 2 {Bt), V^'^iBr), and V^'^iBr) as defined in [19]. 

We also introduce Wp{B), where Z > 0 is not an integer, because initial data will be taken 
from these spaces. The space Wp{B) with nonintegral / , is a Banach space consisting of elements 
of Wp^ ([Z] is the largest integer less than Z) with the finite norm 


llp.R — (u)p g + 


(0 

p,B 


l|(W) 

\\p,B 


where 


Mp!B=l^\\d>\\p,B 
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s=0 


and 



dyu{y)\^ 


dy 

y\n+p{l-ll]) 


_1 

P 


I I 

Wp^{dBT) spaces with non integral Z also play an important role in the study of boundary 
value problems with nonhomogeneous boundary conditions, especially in the proof of exact 
estimates of their solutions. It is a Banach space when p > 1, which is defined by means of 
parametrization of the surface dB. For a rigorous treatment of these spaces, we refer the reader 
to page 81 of Chapter 2 of [19]. 

The use of the spaces Wp ^ (dBr) is connected to the fact that the differential properties 
of the boundary values of functions from Wp’^{Bt) and of certain of its derivatives, d^dj, can 

I i. .. 

be exactly described in terms of the spaces Wp^ {BBt), where 1 = 2 — 2r — s — K 

For 0 < a,/3 < 1, C°‘’^{Bt) is the Banach space of Holder continuous functions u with the 
finite norm 

{x,t)^J3T 


where 


u 


(a)_ 

x,13t 


sup 

{x,t),(x' ,t)£l3T 

x^x' 


\u{x,t) — ^(a;',^)] 
ja; — a;'j“ 


and 




sup 

{x,t),{x,t')^BT 


lu(a:, t) — u{x, t')\ 


We shall denote the space C^’^{Bt) by C^{Bt)- is the set of all continuous 

functions u : Bt —>■ K", and is the set of all continuous functions u : Bt ^ 

K” for which Uxf is continuous for all 1 < i < n. C'^’^{Bt,^'^) is the set of all continuous 
functions u : Bt —>■ R" having continuous derivatives UxiTUx^xj and Ut in Bt- Note that similar 
definitions can be given on Bt- Moreover notations and definitions for Holder and Sobolev 
Spaces on manifolds are similar to the ones used in the Handbook of Global analysis [17]. More 
developments on Sobolev spaces, Sobolev inequalities, and the notion of best constants may be 
found in [6], [7], [12] and [34]. 
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2.2. Preliminary Estimates. For completeness of our arguments, we state the following 
results, which will help us obtain a priori estimates for the Cauchy problem on the manifold 
M, and prove the existence of solutions in Lemmas 2.1, 2.4 and 2.6 can be found 

on page 341, Chapter 4 in [19], as (2.24) and (2.25) on page 49 in [18], and [12] respectively. 
Lemma 2.2 is stated as Lemma 3.3 in Chapter 2 of [19]. 

Let B he a bounded domain in R™ with smooth boundary dB belonging to the class 
with /i > 0 such that B lies locally on one side of the boundary dB. Let T > 0 and p > 1. 
Suppose 0 G Lp{Bt), Wo G Wp{B), 7 G Wp^~p'^~^{dBT). Also, let the coefficient matrix 
(oij) be symmetric and continuous on Sr, and satisfy the uniform ellipticity condition. That 
is for some A > 0 


aij{x,t)^i^j > Aj^j^ for all {x,t) G Bt and for all ^ G 

* j'=i 


Finally, let the coefficients be continuous on Bt- Consider the problem 


dw ^ . , d'^w 


dw 


ij'=l 


dxidxj 

i—1 


+ '^ai{x,t)—= Q{x,t) {x,t) e Bt 


( 2 . 1 ) 


dx_ 

w = 7 (x, t) {x, t) G dBT 
'^\t=o ='^o{x) x€B 


Lemma 2.1. Let p > 1 with p ^ and in the case p > |, assume the compatibility 
condition of zero order, wo\dB = 7|t=o- Then (2.1) has a unique solution w G Wp^'‘^{Bt), and 
there exists C > 0 depending on T,p and B, and independent ofQ,wo and 7 such that 

IkllS. ^ ^(l|0|lr,B. + ll^ollp^+ \h\\p,ot~^^^) 


Lemma 2.2. Suppose q > p, 2 — 2r — s — {m + 2) > 0 and 0 < S < min{(i; Vt}. 

Then there exists ci,C 2 > 0 depending on r,s,m,p and B such that 


\DIDI 




\p,Bt 


for all u G (Sp). Moreover, if2—2r—s— > 0, then for 0<a<2 — 2r—s — 

there exist constants 03,04 depending on r,s,m,p and B such that 




< C 3 S 


||( 2 ) 

llp,0T 


C 4 S 


-(2r+s+h2^+a)| 


Ip.Bt 


for allu€W^’\Br). 


Corollary 2.3. Suppose the conditions of Lemma 2.1 are fulfilled and p > . Then 

there exists c > 0 depending on m,p and B such that the solution of problem (2.1) is Holder 
continuous, and 


(2_ m±2) 


< cllwll 


( 2 ) 

P,Bt 
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Lemma 2.4. Suppose 1 < p < oo. If p < m then Wp{B) embedds continuously into 

Wp ^\dB) and Lq{B) for p < q < p* = 7 ^^- Furthermore, if e > 0 there exists > 0 such 
that 


IKb < 4vX,b 


c.lkll 


p 

1,B 


for all V € Wp[B), and 


Ikll 


2 

2MB 


< e||Wx||2,B 


+ c.|k|| 


2 

2,R 


for all V G Wf{B). 

Lemma 2.5. Let p > m and 0 < a < 1 — ^. Then Wp{B) embedds compactly in C°‘{B). 

Lemma 2.6. Let M be a compact Riemannian manifold of dimension m > 1 and p > m. 
Then the embedding Wp{M) C C°'{M) is compact for all 0 < a < 1 — ^. 

The following result follows from the Gagliardo Nirenberg inequality in [ 8 ] on bounded 
domains, and Young’s inequality on page 40 in [18]. 


Lemma 2.7. Let e > 0 and 1 < p < 00 . Then there exists C^^p > 0 such that 




for all V G Wp{B). 

3. Statements of Main Results. The primary concern of this work is the system 


(3.1) 


ut = DAu + L{{u) 
vt = DAmv + F{u,v) 

U = Uo 


V = Vo 


X 0 <t <T 

X e M, 0 <t <T 

X e M, 0 <t <T 

X € fl, t = 0 
X G M, t = 0 


where D and D are kxk and mx m diagonal matrices with positive diagonal entries, F = (Fi) : 

XK™ K™, G = (Gj) : XK™ and H = (Hj) : ]R'= and uq = (uoj) G W^{n), 

Vo = (voi) G Wp{M) with p > n. Also, uo and vo satisfy the compatibility condition 


= G{uo,vo) 


on M. 


Remark 1 . Since p > n, uo and vq are Holder continuous functions on Q and M respec¬ 
tively (see [1], [(§]). 

Definition 3.1. A function {u,v) is said to be a solution of (3.1) if and only if 
uGG{n X [0,T),R'=)nC^’°(ll X (0,T),M'=)nC^’^(O X (0,r),R'’’) 
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and 

vGC{Mx X (0,T),R'") 

such that (m,t) satisfies (3.1). IfT = oo, the solution is said to be a global solution. Moreover, 
a solution {u,v) defined for 0 < t < 6 is a maximal solution of (3.1) if and only if (u,v) solves 
(3.1) with T = b, and if c? > 5 and (ft, v) solves (3.1) for T = d then there exists 0 < c < b such 
that (M(-,c),n(-,c)) ^ {u{-,c),v{-,c)). 

We say F, G and H are quasipositive if and only if > 0 whenever ^ G and 

( G with = 0 for i = 1, ..., m, and > 0, 77j(C) > 0 whenever ^ G R™ and C G R+ 

with Q = 0, for j = 1,..., k. 

The purpose of this study is to give sufficient conditions guaranteeing that (3.1) has a 
global solution. The following Theorems comprise local and global existence of the solution. 

Theorem 3.2. Suppose F, G and FI are locally Lipschitz. Then there exists T^ax > 0 
such that (3.1) has a unique, maximal solution {u,v) with T = T^ax- Moreover, if < oo 
then 


limsup ||u(-,t)||oo,n + limsup ||u(-, t)||oo.M = cc 

i i ^^max 

In addition to the assumptions stated above, we say condition Vij holds for 1 < j < and 
1 < i < m if and only if 
{Vijl) There exist a, fi, a > 0 such that 

(jFi{C,,v)+Gj{C,,v) < a{C,j+v^+l) and i7j(C) < /3(Ci+l) for all v G R>o, ^ G R|o 
(I^ j 2) There exists Kg > 0 such that 

Gj (C, v) < Kg{Q +v, + \) for all ^ G R^q, C G R|o 
(I^ j3) There exists I G N and Kf > 0 such that 

F^{C, iz) < Kf{\C\ + + 1)' for all i/ G R^q, C G R|o 

Remark 2. (I^j2) is related to the so - called linear “intermediate sums” condition used by 
Morgan in [21], [22] in the special case when the system has only two equations. This condition 
in [21], [22], as well as [23] pertains to interactions between the first m-1 equations in an m 
component system. Again, see [21], [22] and [23]. (14 jT) helps control mass, and allows higher 
order nonlinearities in F, but requires cancellation of high-order positive terms by G. (I4j3) 
implies F is polynomially bounded above. 

Remark 3. We will show that (14^1) provides Li estimates for Uj on 12 and M, and Vi 
on M. (t4j2) helps us bootstrap Lp estimates for Uj on M x (0,Tmax) o,nd 12 x (OjTmax); cLiT'd 
Vi on M X {0,Tmax)- Finally, (14^-2) and (t4j3) allow us to use Lp estimates to obtain sup 
norm estimates on Uj and vt. 

Theorem 3.3. Suppose F, G and IT are locally Lipschitz, quasi positive, and uo,vo are 
componentwise nonnegative functions. Also, assume that for each 1 < j < fc and 1 < i < m, 
there exists h G {l,...,fc} and kj G {l,...,m} so that both t4_;^ and satisfied. Then 

(3.1) has a unique component-wise nonegative global solution. 

Corollary 3.4. Suppose k = m = 1, F, G and H are locally Lipschitz and quasipositive, 
and uoyVo are nonnegative functions. //Vi,i is satisfied, then (3.1) has a unique nonnegative 


V. SHARMA AND J. MORGAN 


global solution. 

In the process of obtaining our results, we will derive Wp’^{Mt) estimates of the Cauchy 
problem on M^, and Holder estimates of the solution to the Neumann problem on Hx. The 
Holder estimates for the solution to the Neumann problem are given as a comment in Brown 
[4]. We give the statement as Theorem 3.6 below, and supply a proof in section 5. Let d,d > 0. 
Consider the systems 



'ft 

— d^M'f + / 

(e,t)GMx(o,r) 

(3.2) 

4'l 

lt=o 

= 4'o 

e G M 

and 





Tt 

= d/S.p + 9 

xGH, 0<t<T 

(3.3) 

dp 

= 1 

X e M, 0 <t <T 



= Po 

X € n, t = 0 


Theorem 3.5. If 1 < p < oo and T > 0, then there exists Cp,T > 0 such that whenever 
2-i- 

'I'o £ Wp ^ (M) and f £ Lp{Mt), there exists a unique solution 4' G Wp'^^Mx) of {3.2), and 


Theorem 3.6. Suppose p > n + 1 and T > 0 and 9 £ Lp{V, x (0, T)), 7 G Lp{M x (0, T)) 
and ipQ G Wp{IX) such that 

d-^^="f{x,0) on M. 
or] 

1 1 

Then there exists Cp^T > 0 independent of 9, 7 and po o-n-d a unique weak solution p € ¥ 2 '^ (^t) 
of (3.3), such that if 0 < /3 < 1 — then 

^ Cp^mka^ + wiWpMr + ii^oiiS) 


The proofs of Theorems 3.5 and 3.6 are given in sections 4 and 5. The remaining results are 
proved in section 6, and examples are given in section 7. 

4. estimates for the Cauchy problem on a manifold. Let n > 2 and M be 

a compact n — 1 dimensional Riemannian manifold without boundary. Consider (3.2) where 

2-^ 

d > 0, / G Lp{Mt) and H/q G Wp ^ {M). Searching the literature, we surprisingly could not 
find Wp'^{M t) estimates for the solutions to (3.2). Tracing through the work in this direction, 
we found that Huisken and Polden [15] and [25], and J.J Sharpies [31] give a result in the 
setting where p = 2. Using their estimate, we obtain Wp^^{Mt) a priori estimates 

for solutions of (3.2) for all p > 1. For a > 0 and smooth functions f,g : M x [0, 00 ) — >■ K, 
Polden considered weighted inner products: 



e ^“‘(/(•,t),5(-,t))i2(M)dt 


if, 9 )ll. 
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if, 9) 


LW^ — 




{fi9)LW^ — e 

Jo 

{f,g)wWa = + {Dtf,Dtg)LLa 

Where LLa, LWa and WWa are the Hilbert spaces formed by the completion of C°°{M x [0, 00)) 
in the corresponding norms, and WW^ is the completion of subspace of C°°{M x [0, 00)) with 
compact support in WWa- See [31] for the proof of the following result. 

Theorem 4.1. Suppose 4*0 lies in W^iM) and f € LLa{M x [0,oo)). Then for sufficiently 
large a, the system (3.2) has a unique weak solution in WW^- 


Furthermore using a priori estimates in [31], they showed that the solution belongs to W^’^iM x 
[0,oo)). 

Theorem 4.2. Let T C WWa be the unique solution of (3.2) with Ibo G Wf{M) and 
f € LLaiMx)- Then T £ LVVf, and there exists C > 0 independent o/Tq and f such that 

miw^<cmorwi^M) + \\f\\iLj 

Proof. See Lemma 4.3 in [31]. □ 

The result below is an immediate consequence. 

Corollary 4.3. Let 0 < T < 00 . Suppose Tq £ Wf{M) and f £ L 2 {Mt)- Then there 
exists a unique weak solution to (3.2) in W 2 '^{Mt), and there exists C > 0 independent o/Tq 
and f such that 

+ ll/lli2(MT)) 

We will use the W 2 '^{Mt) result to derive W^d{MT) a priori estimates for solutions to 
(3.2) for all p > 1. To obtain these estimates, we transform the Cauchy problem defined locally 
on M to a bounded domain on and obtain the estimates over this bounded domain. 

Then we pull the resulting estimates back to the manifold. Repeating this process over every 
neighborhood on the manifold, and using compactness of the manifold, we get estimates over 
the entire manifold. 

Let be a subset of K+ with following property: 

2-^ 

p > 1 belongs to P if and only if there exists Cp^r > 0 such that whenever Ibo £ Wp ’’ (M) 
and / £ Lp{Mt), then there exists a unique T £ Wp’^^Mx), such that T solves (3.2) and 

Note: From Corollary 4.3, 2 G P. Also note that we can prove Theorem 3.5 by showing 
P = (1,00). 
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Lemma 4.4. [2,oo)cJ'. 

Proof. We will show that ii p G 7' then [p,p + C 7. To this end, let p G 7 

2 -- 

and q G [p,p+ such that Tg G Wq “ {M) and / G Lq(Mr). Then / G Lp^Mx) and 
2 -.^ 

4^0 S Wp ^ (M). Since p G 7, there exists Cp^x > 0 independent of d'g and /, and a unique 
4^ G Wp’^{Mx) solving (3.2) such that 

( 4 - 1 ) < Cp,xmpMr + 

Let B{0, 1) be the open ball in of radius 1 centered at the origin. Now, M is a 

manifold. Therefore, for each point ^ G M there exists an open set Vj of M containing ^ and 
a diffeomorphism cj)^ : B{0, 1) Vj. Let 4> = 41 o ((ij, / = f o and $o = 4'o o cj)^. Using 
the Laplace Beltrami operator (defined in [30]), (3.2) takes the form 

yj det g d,<^) + f{x, t) xG 5(0,1), 0 <t<r 

Vdet g 

(4.2) 4> = 4>o xG 5(0,1), t = 0 

where g is the metric on M and is the i, entry of the inverse of the matrix corresponding 
to metric g. That is, in the bounded region 5(0,1) x (0,T), we have 

n—1 n—1 

(4.3) £($) = ^ + X! = / 

i,j—l i—1 

(4.4) = $0 
where. 


Gij = d g*-’ 


«i = rj^ dj{g''Wdet g) 
y/det g 

Note 4/ G Wp'^{Mx) implies 4> G Wp'^{B{Q, 1) x (0,T)). Take 0 < 2r < 1 and define a cut off 
function if G [0,1]) such that. 


(4.5) 


V’(x) 


1 Vx G 5(0, r) 

0 Vx G R"-i\5(0,2r) 


In Q = 5(0,2?’), Qx = 5(0, 2r) x (0,5) and Sx = 95(0, r) x (0,5), w = satisfies the 
equation 


dw 

'm 


n—1 

i j'=l 


d'^w 

dxidxj 


+ E 



w = 0 




(x,t) G Qx 

(x,t) G Sx 
t = 0,x G Q 
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where, 


n—1 r, , n—l 

0 = /V'- 2 ^a,, —^ 


ij=l 


ay 

dxjdxi 




Since V' G [0,1]) and $ G M^p’^(i?(0,1) x (0,T)), therefore 6 — f'ip £ Wp’^{Qt)- 

Case 1. Suppose p < n. From Lemma 2.4, 6 — ftp £ ^niin{q p+-Ed_}(*5T)- hr particular 

1 2 - ^ 
since p + < p + and fip G Lq^Qx), we have 6 G LqiQr)- As a result 

II^IU.Qt < II/V’II^.Qt + C'i||4)||q,_Qj, + C'2||$3 ;||,_Qj, 

< \\mq,QT+cmu,Q.+c2\m\%^ 

where Ci,C 2 >0 are independent of /. Now in order to estimate ||‘i>a,||p^Q^, apply the change 
of variable 




( 1 ) 

P,Qt 


||4/,|det((0^-i)')||| 


( 1 ) 

p.(0«(Q)) 


T 


and using (4.1), we get 

\\0\\q,Q^ < mWq.QT + CmkQ. + C2p,T{\\f\\pMT + Il'I^ollpV^) 

where C 2 p^T > 0 is independent of / and Iho. At this point, we need an estimate on ||4>||g^Q,j,. 
Again ||$||q.Qj, = ll^'l det((()ij"^)')|||q,(,^^(Q))^ and from Lemma 2.4, 

||vI/|det((<(.r')')lll,.(,,(Q)), < C|l'I^|det((<(.r')')lllS<^,(Q))^ 


Thus 

(4.6) \\e\\q,Q^ < Kp^TiWfWpMT + 114'ollpV^) 
where Kp^x > 0 is independent of / and Tg. 

Since gij are functions on the compact manifold M, aij and Ui satisfy the hypothesis 
(bonnded continuous function in Qx) of Lemma 2.1. Therefore using Lemma 2.1, 

(4.7) lkllS,<t^,.T(||0|U.Q. + ||V^4>o||^^-"^) 

where Cq^x > 0 is independent of 9 and ip^Q. Combining (4.6) and (4.7) we get, 

< KpAWfWpMT + + IIV'^oII^'q^’) 

where Kp^x > 0 is independent of /, 9 and ^$o. Note that w = $ on Wx = B{0,r) x (0,T). 
Thus 

< Kp,x{\\f\\pM. + + mo\fq,Q~^^) 


(4.8) 
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Observe (4.8) is over B{Q,r) x (0,r) C M” ^ x K+. To get the estimate back on the manifold, 
apply the change of variable, = ||'I'| det(((/)“^) )| 11 ^^ 0 ( 14 /.^,) and using first mean value 

theorem of integration there exist ^ G ^(Wr), and j, ^ such that 

(4-9) l|4^lliVr) ^ + l|4^o|ll';|)) 

So far, an estimate in one open neighborhood of some point ^ G M is obtained. As 
one varies the point ^ on M, there exist corresponding open neighborhoods Vj and a smooth 
diffemorphisms cj)^ : 5(0,r)—which results in different for every Vj. Consider an 

open cover of M such that M = U{eM 4^- Since M is compact, there exists {^i, ^ 2 , ^n} such 

that M C\J V'jj and K j .corrresponding to each . Let, Cp^M,T = J2i<i<N ^p t ( ■ 

1<3<N 

Inequality (4.9) implies 

\m[%<CpMM\\f\u,M. + \\'i'oCM^^) 

Thus [p,p+ C 5. 

Case 2. Suppose p > n. 

2 -^ 

By Lemma 2.4 and Theorem 4.12 in [1], if g G [p, 00 ), dlo G Wq " {M), and / G Lq{MT) 
then 6 G Lg^Qx), and proceeding similarly to Case 1, we get 

ll'I^lli'k<C.,M,T(||/||,.M, + ||4/o||jy^) 

where Cq^M,T > 0 is independent of /, 9 and Hence [2, 00 ) C TH 

Proof of Theorem 3.5: From Lemma 4.4, we have [2, 00 ) C T. It remains to show that 

(1,2) C J'. Let 1 < p < 2 , / G Lp{Mt) and 4'o G W^~i{M). Since C°°(M^) is 

_ 2-1 

dense in Lpi^Mx) and C°^{M) is dense in Wp ’’(Af), there exist a sequences of functions 
{/fc} C C°°{Mt) and {'I'ofe} C C°°{M) such that fk converges to / in Lp{Mt) and 4'ofc 
2-2 

converges to 4>o in Wp ” (M). Define a sequence {4'*,} such that, 

Tfc, = + /fc ^ G M, 0 < t < T 

(4.10) 4'fc = ^ofc /eM, t = 0 

Now, transform system (4.10) over a bounded region in Similar to the proof of Lemma 

4.4, for each point ^ G M there exists an open set oi M containing ^ and a diffeomorphism 
(j)^ : 5(0,1) Vj. Corresponding to each k, let fk = fk ° 4’ofc = 'I'ofc o and using the 
Laplace Beltrami operator, (4.10) on 5(0,1) C 5 takes the form 

(4.11) = + xG 5(0,1), 0<t<T 

y/det g 

$fc = 4>0fc xG5(0, 1 ), t = 0 

Consequently, in a bounded region 5(0,1) x (0,T) of the Euclidean space, we consider 

(4.11) in the nondivergence form defined in (4.3) for each $fc, with / replaced by fk and $0 by 
‘J’Ofc- 
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Taking 0 < 2r < 1, using a cut off function ip € [0,1]) defined in (4.5), and 

defining Q = B{0,2r), Qt = B{0,2r) x (0,T), St = dB{0,r) x (0,T), and Wk = ip^k, we see 
that 


where, 


dwk 

~dr 


n—1 

i j'=l 


d'^Wk 

dxidxj 


+ E 


dwk 
' dxi 


Ok 


Wk = 0 

^fc|t=0 = V'^Ofe 


{x,t) G Qt 
{ x,t) G St 

t = 0,x € Q 


n—1 

Ok = fk'>P - 2 ^ ( 

2=1 


d^k dip 

dxi dxj 


n — 1 

- ^ a* 

ij'=l 


dxidxj 


a. 

2 = 1 


Note that fk and 4'ofc are smooth functions. Therefore Lemma 4.4 guarantees G Wg’^{QT) 
for all q > 2. Thus Ok G Lq{QT) for all q > 2. Recall ip G [0,1]). Using Lemma 2.7 

for e > 0, there exists > 0 such that 

||^fc||p,QT ^ ll/fc^llp,QT + -^lll^fcllp.Qr + ^2\\^kx\\p,QT 
— ll/fe|lp,QT ff -^lll^fcllp.QT 

(4.12) + M2{€\\^kxx\\p,QT + ^(W^kWp^Qr) 

Here Mi, M 2 > 0 are independent of / and d'o- At this point we need an estimate for ||$fe||p^Q,j,. 
From Lemma 2.4 forl<p<n<(7 there exists Ce > 0 such that 




p 

L pq (Qt) 

q — p 


< e(ll^fea;|| 


p 

P,Qt 


1$ 


kt 


l^,Q. 


) + C.|l4>fc|| 


P 

1 ,Qt 


Since p < from Holder’s inequality, e and get scaled to e > 0 and Cg > 0 (with e —>• O'*' 
as e —>■ 0+), and 


(4.13) 


From (4.12) and (4.13), 


ll^fcllp.Qr ^ ^(ll^fctllp.QT + ll‘^fcx|lp.QT) 


ll^fellp.QT ^ (Afl + M2Ce){e{\\^kt\\p,QT + ll‘J’fca:|lp,QT) + C'gjl || ) 

+ \\fk\\p,QT + M2e\\^kxx\\p,QT 

Recall gij are functions on the compact manifold M. Therefore Oij and at satisfy the 
hypothesis (bounded continuous function in Qt) of Lemma 2.1. Using Lemma 2.1 for p |, 

(4.14) IkfellJE ^ ^pTdl^fcllp.Q. + mok\fp,Q^^) 

where Cp^T is independent of 0 and ip^Q. Combining (4.12) and (4.14), we get 

||w/c||p^Q.j, < Cp^Ti\\0k\\p,QT + 

< Cp,T{\\fk\\p,QT +M2e||4> kxx\\p,QT 

+ {Ml + M2Ce){e{\\^kt\\p,QT + l|4’fexllp,QT) + C'g|l$fc||^ g^) 

+ mokCg^^h 
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Note that Wk = on Wt = B{0,r) x (0,T). Thus 

Cp,T{\\fk\\p,QT + M2e\\^kxx\\p,QT 

+ (Ml + M2Ce){e{\\^kt\\p,QT + 11‘i’fcxllp.QT) + 

(4.15) +mok\\p,Q^^} 

Observe (4.15) is over 5(0,r) x (0,5) C x R_|_. To get an estimate on the manifold, 

apply the change of variable, = ||4'fe| 'let((^“^)')|||^^^(^,^,^^ and using first mean value 

theorem of integration there exist ^ € 4>{Wt), and j, ^ such that 

II4^^lip,^(W t) — ^P>?5{ll/fcllp,<^’(QT) + .^2e||'I'fca:x llp,<;i(QT) 

+ (Ml + M2Ce)(e(||'I'fct||p_,^(Q.j,) + ||'I'fc2,||p_0(Q,j,)) + C'?l|4'fe|li,(0(Qr))) 

(4.16) +l|4^0fe||p';(|)} 

So, an estimate in an open neighborhood of a point ^ € M can be obtained. As one varies the 
point ^ on M, there exist corresponding open neighborhoods and a smooth diffemorphisms 
(j)^ : 5(0,r) —V^, which result in different for every V^. Consider an open cover of 

M such that M = (J^eM 4^- Since M is compact, there exists {■^1,^2, such that M C 

U V^. and C f corrresponding to each Vj^.. Let Cp^r = J2i<i<N^pi t- Inequality 

l<j<N 

(4.16) implies 

I|4^^IIp,Mt — C'p.tIIIA-IIp.Mt +M2e||\E' kxx lIpjAfr 

+ (Ml + M2Ce){e{\\'^ ktWpMT + \\^kx\\p,MT) + 

(4.17) +\\'fok\\p,;/h 
Also, a simple calculation gives 

ll'hfclll.Mr < IIAIIi.Mt + l|4'0fc||l,M 

Now, choose e > 0 such that. 


max{C'p_TM2e, Cp^T^Mi + M2Ce)} < ^ 

For this choice of e, (4.17) gives the estimates 

Il'l'fclli^U < C’pTdl/fclUM. +t^.(IIMIl,M. + ||4/0fc||l.M) + \\'i^0k\&) 

(4.18) l|4^felli'k < kp,Ti\\fk\\p,MT + Il'I'ofcllpV^) 

where Kp^x > 0 is independent of fk and IPofe. H remains to show that the sequence {d/fc} 
converges to a function in Wp’^{Mt), and 4' solves (3.2). From linearity and (4.18), if 
TO, 1 G N then 'I'm — 4// satisfies 

(4'm-4'z)4 =jAM(4'm- 
4'm - 4-/ = 4-0^ - 


4'i) + fm — fl 


c G M, 0 < t < 5 
C G M, t = 0 
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and 

Il'I'm ~ ^ Kp^TiWfm ~ fl\\p,MT + Il'^'Orn ~ ) 

This implies {'I'fe} is a Cauchy sequence in Wp’^^PIr), so there is a function ^ G Wp’^{Mt) 

2-^ 

such that —>• vp. Then fk converges to / in Lp^Mr), 4'ofc converges to 4'o in Wp ^(M), 

and converges to 'h G Wp’^{Mt)- Therefore 'I' solves (3.2), and (4.18) implies 

m%^<KpM\\f\\p,MT + \\'i'o\\p,J^) 

Hence T = (l,oo), and the proof of Theorem 3.5 is complete. 

5. Holder Estimates for the Neumann problem. The following result is a version of 
Theorem 9.1 with Neumann boundary conditions, referred to in chapter 4 of [19] on page 351. 

Lemma 5.1. Tetp > 1. G Tp(r2 X (0,T)), (/3o G IHp ^\^)andjGWp ^’’(Mx 

(0,T)) with p 3 . In addition, when p > 3 assume 



on M X {0} 


Then (3.3) has a unique solution (p G Wp'^(yi x (0,T)) and there exists C dependent upon 
VL,p,T, and independent of 0,ipo and 7 such that 


ll‘^llp.(nx(o.T)) — ^(ll^llp.(nx(o.T)) + llv^ollp^o” + 


(2-f) 


fi-i i-J-i 

p’2 2p> \ 

llp,(9nx(o,r))l 


if 


Definition 5.2. p is said to be a weak solution of system (3.3) from V 2 (Hy) if and only 


r f 

lo Jq Jo JdQ iJv 

= f iyix,0)ip{x,0) 

Jq 


d \7ix.\7(p — 


Ov 


/O JQ 


lo JQ 


for any v G that is equal to zero for t = T. 


We also need a notion of solution of (1.3) which was first introduced in the study of Dirichlet 
and Neumann problems for the Laplace operator in a bounded domain by Fabes, Jodeit 
and Rivier [9] . They used Calderon’s result in [5] on continuity of Cauchy integral operators 
for curves. Further in [10], Fabes and Riviere constructed solutions to the initial Neumann 
problem for the heat equation satisfying the zero initial condition in the form of a single layer 
heat potential, when densities belong to Lp{M x (0,r)), 1 < p < 00. We will consider the 
solution to (1.3) in the sense of one which is constructed in [10]. 

The following result plays a crucial role for that construction of solution to make sense, 
and is proved in [10]. 


Proposition 5.3. Assume Q, is a domain and for Q G M, riq being the unit outward 
normal to M at Q. For 0 < e < t set 


Mf){Q,t)= f 
Jo 


t M 


{y-Q,VQ) 

{t — s)^+^ 


exp - 


\Q-y\^ 

4(t- s) 


/(s, y) da ds 


Then 
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1. For every 1 < p < oo there exists Cp > 0 such that \Je{f){Q,t)\ = J[f)(Q,t) 

satisfies 

\\J{f)\\Lp{Mx{o,T)) < C'pII/||lp(Mx(o,t)) for all f € Lp{M x (0,r)) 

2. linif .(1+ JJf) = Jif) exists in LJM x (O.T)) and pointwise for almost every (Q.t) € 
(M ^(0, T)) provZ f e L,(M X (0, T)); 1 < p < 

3. CnI + J is invertible on Lp{M x (0, T)) for each 1 < p < oo and Cn 0. 

We consider the case d = 1 below. The extension to arbitrary d > 0 is straightforward. For 
Q £ M, {x, t) £ FIt and t > s, consider 

exp 

W{t - s,x,Q) = - ; ^ and g{Q,t) = -2[-c„/ + J]~^j{Q,t) 

{t-s )2 

where Cn is given in [10]. 

Definition 5.4. tp is said to be a classical solution of system (1.3) with d = 1 and, 
7 £ Lp{M X (0, T)) for p > 1 if and only if 


p{x,t) 



s, X, Q)g{Q, s) da ds for all (x, f) £ 


Remark 4. When 0 = 0 and ip{x,0) = 0, the weak solution of (3.3) is the same as the 
classical solution of {1.3). 


In order to prove the classical solution p to (1.3) is Holder continuous, let {x,T), {y,T) £ 
such that 

p{x,T) = ( f W{T — s,x,Q)g{Q,s) da ds 

Jo Jm 


and 


F{y,T)= / W{t - s,y,Q)g{Q,s) da ds 

Jo JM 

Without loss of generality we assume 0 < r <T. Consider the difference 

tp{x,T)-ip{y,T)= f f {W{T - s,x,Q) -W{t - s,y,Q))g{Q,s) da ds 
Jo Jm 

+ f [ W{T — s,x,Q)g{Q,s) da ds 
Jt Jm 


Lemmas 5.5, 5.6 and 5.7 provide estimates needed to prove p is Holder continuous. Throughout 
the proofs we assume p' = . 


Lemma 5.5. Letp > n+1. Suppose {x, T), {y, r) £ Ut with 0 < r < T and TZ‘^ = {(Q, s) £ 
M X (0,t) : |x — QI + jT — sj^ < 2(|a; — y| + |T — r|^)}. Then for 0 < a < 1 — there exists 
Ki > 0 depending onp,n,fl, T and independent of g G Lp{M x (0,r)) such that 

f \{W{T - s,x,Q) -W{t - s,y,Q))g{Q,s)\ da ds 
Jn‘= 


<i£i (^|a;-?/| + |T-r|2^ II 5 llp,Mx[o.r] 
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Proof. 


j \{W{T - s,x,Q) - W{t - s,y,Q))g{Q,s)\ da ds 

exp ( ~4'(T-y ) exp 



\g{Q, s)| da ds 

(r-s)': 



Ilffllp.TC' 


By hypothesis p > n + 1. Pick 0 < e < ®et N = ” ^ Then there exists c > 0 such 


that w" ■ exp(—w) < c • for all re > 0. Consequently, 


exD f 

«xp 4(T-s) 

/tc<= (T-s)^ 


2\ \ ZT 


exp 


-p'\y-Q\ ^ ^ ^ 

4(r-s) 




< 


< 


■■N 


'n‘^ (r —s)’ 
^ / 


\\9\\p,n<= 




N 


c-N 


Cl ( / ds 

0 


1 


’■R‘ 


f 1 


— da 


N 


\\9\\p.n<= 




da 


9 llp,7?.= 


where A = {Q G M : \x — Q\ < 2\x — y\ + \T — t\^}. Since |T — t| < \T — s\,TZ^ C A x (0,r). 
Let py = \y — Q\, px = \x — Q\. Notice that in A, 0 < < ‘2\x — y\ + \T — rjs and 

0 < py < |a; - 2/1 + Pa, < 3|a; - j/| + |T - r|^. Therefore, 


II, . n — l — e — np' 

Cl [ {t — s) 2 ds 


II, . Ti — l — e — np' 

+C 2 I / (T-s)-5- ds 


1 


da 


< 


Cl 


n-l-e-np' 

[T — Sj 2 ds 


L \y-Q\^-^-^ 
|x-Q|— 1 - 

p3|x-y| + |T-r|i 


da 


9 lip,7?.= 


+C2 I 1^ (T-s) 


X — 1 — e — np' 


f2\x-y\ + \T-T\ 2 


ds 


r^-^dr 


^dr 


IIsIIp.tc^ 
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a^ . n + l-e-np' 


—M 

^3|a;-?/| + |r-r|2j 


C‘2 f n-\-\ — e — n-p' _ . n + l — e — np' \ "p" / , , , , 1 \ 

+^{t = -(T-r) 2 j (2\x - y\ + \T - T\-^j 


Il5llp.7^= 


By hypothesis, p' < . Therefore, there exists Ki > 0 depends on p,n and T such that 



s,x,Q) - W{t - s,y,Q))g{Q,s)\ da ds 


< Ki [\x - y\ + \T - ” II 5 IIp.mx[o.t] 


The result follows since 0 < e < is arbitrary.D 

The proof of the following Lemma makes use of Brown’s corollary to Theorem 3.1 in [4]. This 
also provides a proof for the remark made in [4] after Lemma 3.4. 

Lemma 5.6. Let p > n + 1. Suppose (a:, T),{y, r) £ LIt and TZ = {{Q, s) G M x (0, r) : 
2{\x — y\ + \T — t\^) < |x — (51 + |T — s|5}. Then for 0 < a < 1 — there exists K 2 > 0 
depending onp,n,Q, T and independent of g £ Lp(M x (0,T)) such that, 


f \{W{T - s,x,Q) - W{t- s,y,Q))g{Q,s)\ da ds 
Jn 

< K 2 (|a; -y\ + \T - t|^) || g ||p,Mx[o.r] ■ 


Proof. Using the Theorem 3.1 in [4], we have 


[ li^iT - s,x,Q) - W{t - s,y,Q))g{Q,s)\ da ds 
Jn 

< f C (——(1 + (T- s)^)exp ( ) |g((5,s)| da ds 


In \ |r - s| 2 + |a; - (51 


4(T-s) 


<^1(2 


<^i 


m 


5k y |T - s |2 + |a; - (51 / \T-s)^ 

1 exp[±^] 

:-0|“ (T-»)T \<HQ,s)\d^d, 


where Hi = (7(r2 + 1) and i)i = Hi (i)^ “ ^|T — t| 2 + |a; — ?/|^ . By hypothesis, n + 1 


(n + a)p' > 0. Pick 0 < e < (n + 1) — (n + a)p' and set N = 


n—l—e—ap 


. Then there exists c > 0 
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such that zu™ • exp(—zc) < c - N for all zc > 0. Consequently, 


D, 


<Di 


1 exp 


V 4(t-s) J 


I 


da ds \\ g ||p,K 


c-N 


’n 




N 


9 \\p,n 


< cDi 


{T-s) 




lo Jm (T-s)^ \x-Q\ 


jrr^dads] || 5(s, Q) IIp,mx[o,t] 


If, . — 1 —e — ap^ — np' 

< cDi [ / (T — s) 2 ds 


I 

J M 


'0 Jm k ~ Q\‘' 

Then by change of variable, there exists C, a > 0 such that 


— l — e ) II 5 IIp,Mx[ 0 ,t] 


ds. 

< CDi (. r 


1 


Im \x - Q\ 


n—l — e ) II d llp,Mx[0,r] 


i—^dr 1 \\ g ||p,mx[o,t] 


The result follows. □ 

Lemma 5.7. Let p > n + 1, and suppose (x, T),(z/, r) G LIt- Then for 0 < a < ^ — 
there exists > 0, depending on p,n,Ll and T, and independent of g € Lp{M x (0,r)) such 
that, 



s,x,Q)g{Q,s)\ da ds < K.s{T - tY || g ||p,Mx[r.T] 


Proof By hypothesis p > n + 1. Pick 0<e<n + l — np' and set N = 
exists c > 0 such that ■ exp(—zc) < c ■ N for all w >0. Consequently, 


Then there 


| 1 P(T — s, x, Q)g{Q, s)| da ds 


Jm 

fT p exp 


< 


f -\x-Q\^ \ 

^ ^ -J—^\g{Q,s)\ da ds 


<C3 


T JM [T ~ 


IM 


(T-s) 2 \x-Q\ 


^ 3 Y 3 ^|p(( 5 ,s)| da ds 


<C 3 


{T-Sy '2 II 5 ||p,Mx[r.T] 
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Similarly, by change of variable there exist Ca, a > 0 such that 

< (^3 (^I(T-t)" ' +^| ■ II IIp.Mx[t.T] 

< KsiT — t) 2^^ II g ||p,mx[t.t] 

where K3 > 0, depends on p,n,fl and T, and independent of g € Lp{M x (0,T)). The result 
follows since 0 <e<n + l — np' is arbitrary, and = i — . □ 

Proposition 5.8. Suppose 7 G Lp{M x (0,T)) for p > n + 1. Then the classical solution 
of (1.3) is Holder continuous on Tl x (0,T) with Holder exponent 0 < a < 1 — and there 
exists Kp > 0, depending onp,n,Q and T, and independent of j such that 

\p{x,T) - p{y,T)\ < kp{\T-T\^ + |a; - ?/|) || 7 ||p,Mx(o.T) 

for all (x,T), (y,T) G fix- 

Proof. We prove this proposition for d = 1. The extension to arbitrary d > 0 follows from 
a simple change of variables. Let O be an open subset of 12 with smooth boundary such that 
the closure of 12 is contained in 12. It is straightforward matter to apply cut-off functions and 
Theorem 9.1 in [19] to obtain an estimate for (p in ITp’^(12 x (0,T)). Moreover, there exists 
LpQT independent of 7 such that 


11 ‘^IIp.At ~ '^p.n.rll^llp.MT 

Since p > n + 1, Wpd[fl x (0,T)) embeds continuously into the space of Holder continuous 
functions (see [19]). As a result we have Holder continuity of the solution to (1.3) away from 
Mt- We want to extend this behavior to points near Mt- 

Pick points {x,T), {y,T) G 12t. We know from Fabes and Riviere [10] that the solution of 
(1.3) is given by 


(p{x,T) = [ f W{T — s,x,Q)g{Q,s) da ds 
Jo Jm 


where W{T — s, x, Q) 


f -|x-Q|2 

d 4(T-o) 


{T-s)2 


, g{Q,t) = [I + J] ^jiQ.t) and 


J{g){Q,t) 


lim r 

e^0+ Jo 



{y - Q,vq) 

{t — s)^+^ 



4(2 -s) J 


g{s,y) da ds 


for almost every Q G M (for smooth manifold it is true for all Q), gq being the unit outward 
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normal to M aX Q. 

\ip{x,T) - Lp{y,T)\ = \ j f {W{T - s,x,Q) -W{t - s,y,Q))g{Q,s) da ds 
Jo Jm 

+ [ [ W{T — s,x,Q)g{Q,s) da ds\ 

Jt Jm 

< \ f {W{T - s,x,Q) -W{t - s, y, Q))giQ, s) da ds| 

+ 1 [ iWiT - s,x,Q) - W{t - s,y,Q))g{Q,s) da ds\ 

Jn 

^ Ir (T- s)^)exp -^s) ) 

Where TZ and TlX are given in Lemmas 5.6 and 5.7. Now using Lemma 5.5, Lemmas 5.6 and 
5.7 for 0 < a < 1 — there exists Ki,K 2 , > 0 depending on p, n, fi, T and independent 
of g G Lp{M X (0,T)), such that 


So, 


□ 


\(p{x,T) - (p{y,T)\ < Ki (^|a; - y\ + |T- r|^ j || g ||p,Mx(o.r) 

+ iC 2 (|T-T|^ + II 5 llp.Mx(o.r) 

+ Ko{T - t) 2 ^" II g ||p,mx(t,t) 

\p{x,T) - <p{y,T)\ < kp[\T - t|5 + |a; - ?/|) || g ||p,mx(o.t) 


Now we combine Holder estimates and Theorem 9.1 in chapter 4 of [19] to get the existence 
of a Holder continuous solution to system (3.3) for any finite time T > 0. 

Proof of Theorem 3.6: Chapter 4, Theorem 5.1 in [19] implies (3.3) has the unique weak 
solution. In order to get Holder estimates, we break (3.3) into two sub systems. To this end, 
consider 


(5.1) 


Lp2t = d/^if2 + d 




^>1 = '•Po 


X Gil, 0 <t <T 
X G M, 0 < t < T 
a: G H, t = 0 


ipif. = dAipi 



= 1 



ipi =Q 


a: G H, 0 < t < T 
X G M, Q <t <T 
a; G H, t = 0 


(5.2) 
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From Lemma 5.1 there exists a unique solution of (5.1) in x (0,T)), and a constant 

Cl (T, p) > 0 independent of 9 and i^o such that 


II^2||Sx(0.t) < CiiT,p)m\p,n.io,T) + ll^ll 


djpo 

dr] 


n-i i-U-i 

D ’ 2 2 tj / 


p,(dnx(o,T)) 


) + ll¥>o|lS 


Using proposition 5.8, there exists C' 2 (r, 0) > 0 independent of 7 and po so that the unique 
weak solution to (5.2) satisfies, 


\pii%,T)<C2{T,p) 


llp,Mx(0,r) + l|-^||p,Mx(0,T) 


where 0 < /3 < 1 — By linearity, p = pi + P 2 solves (3.3). Moreover, for p > n + 1, 

W^’^(U X (0, T)) embeds continuously into i^r)- So, there exists C{T,p) > 0 independent 
of 0, 7 and pq such that 

(5-3) l‘^lffi(o,T) - C'(’^iP)(ll^llp,Ox(o,T) + ll7llp,Mx(o,T) + llv^ollp^n) 


Remark 5. lUe will use these Holder estimates to obtain sup norm estimates, and local 
existence results for (3.1). 


6. Proof of Theorems 3.2 and 3.3. 

6.1. Local Existence. 

Theorem 6.1. Suppose F, G and H are Lipschitz. Then (3.1) has a unique global solution. 

Proof Let T > 0, Fix {uo,vo) G Wp{D,) x Wp{Tt) such that they satisfy the compatibility 
condition 

( 6 . 1 ) D^^ = G{uo,vq) on M. 

or] 

Set 


X = {{u,v) G C{n X [0,r]) X C{M X [0,T]) : u{x,0) = 0,V a; G U,?;(x,0) = 0,V a; G M} 
Note (X, II • ||oo) is a Banach space. Let {u,v) G X. Now consider 


Ut = DAU + H{u + uo) 

a: G D, 

0<t<T 

Vt = DAmV + F{u + uo,v + Vo) 

X G M, 

0<t<T 

dU 

( 6 . 2 ) + Uo, + Wo) 

07] 

X G M, 

0<t <T 

U = uo 

X G it, 

t = 0 

0 

II 

X G M, 

t = 0 

From Theorems 3.5 and 3.6, (6.2) possesses a unique 

weak solution 

iu,v) G X 


Wp’^{Mt). Furthermore, from embeddings, {U,V) G C{fl x [0,T]) x C{M x [0,T]). Define 

S : X ^ X via S{u, v) = {U — uq, V — vq), 
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where {U, V) solves (6.2). We will see that S is continuous and compact. Let {u, v), {u, v) G X. 
Then 


Siu, v) - S{u, v) = (U-U,V-V) 
Using linearity, {U — U,V — V) solves 


Ut-Ut = DA{U -U) + H[u + Mo) - H{u + uo) 

X £ n, 

0<t<T 

Vt-Vt = DAm{V - U) + F{u + uo,u + Uo) - F{u + mo,u + uo) 

X £ M, 

0<t<T 

07] 

X £ M, 

0<t<T 

U-U = 0 

X £ 

t = 0 

V -V = 0 

X £ M, 

t = 0 


From Theorem 3.6, if p > n + 1 there exists K independent of H, G, F, u, v, u, v such that 

\\U - Ulloo.nT + 11^ - V\\oo,Mt < K (||U(u + lto,T + Vq) - F{u + Uo,V + Vo)\\p,Mt 

+ \\G{u + Uo,V + Vo) - G{U + Uo,U + Uo)||p,Mt 

+ \\H{u + uo) - H{u + ito)||p,nT) 

Using the boundedness of U and M, there exists K > 0 such that 

\\u - UIIoo.Ot + 11^- V\\oo,Mt < K (||F(u + Uo,V + Uo) - F{u + Uo,V + Vo)\\oo,Mt 

+ ||G(u + uo,v + Vo) - G{u + Moiii + wo)||oo.Mt 
+ \\H{u + uo) - H(u + ito)||oo.nT) 


Since, F, G, H are Lipschitz functions there exists M > 0 such that 

\\U - U||oo,nT + 11^ - V\\oo,Mt < M{\\u - + ||u - w)||oo,Mt) 

Therefore S is continuous with respect to the sup norm. Moreover, for p > n +1, from Theorem 
3.5, 3.6, and Lemma 2.6, there exists C(T,p) > 0, independent of F{u + uo,v + vo), G{u + uo,v + 
vo),H{u + uo), uq and vo such that for all 0 < a < 1 — 0 < /3 < 1 — 

(6-3) \U\al + |U|^], < G(T,p)(||iJ(u + uo)||p,a,r + ||G(u + uo,u + uo)||p,Mr 

+ ||F(u + uo,v + vo)\\p,Mt + Ikollp^L + llwollp^n) 

Using (6.3), S maps bounded sets in X to precompact sets, and hence S is compact with 
respect to the sup norm. Now we show S has a fixed point. To this end, we show that the set 
A={(it,u) £ X : {u,v) = XS{u,v) for some 0 < A < 1} is bounded in X with respect to the 

sup norm. Let {u,v) £ A. Then there exists 0 < A < 1 such that (j, j) = S(u,v). Therefore if 

(u, v) = (u + Auo, V + Avq) then 


Ut = DAu + AF[ {u + Uo) 
vt = DAmv + AF{u + uo,v + Vo) 
du 

F— = AG{u + uo,v + Vo) 

u = Auo 
V = Avo 


X £ fl, 0 < t < T 
X £ M, 0 <t <T 

X £ M, 0 <t <T 

X £ t = 0 
X £ M, t = 0 
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From Theorem 3.6 and H, F and G being Lipschitz, there exists N > 0 such that IKu, '0)11 oo < N, 
with N independent of X,u and v. Since ||('«,n)||oo < ||(m,'0)||oo < N, hence boundedness of 
the set is accomplished. Thus, applying Schaefer’s theorem (see [8]), we conclude S has a fixed 
point {U,V). Further, {U uo,V + vo) is a solution of (3.1). Moreover, bootstrapping the 
regularity of this solution using well known estimates, we obtained a solution to (3.1) according 
to Definition 3.1. 

Finally, we show the solution of (3.1) is unique. Suppose {u,v), {u,v) solve (3.1). Then, 
{u — u,v — v) satisfies 


Ut — Ut = DA{u — m ) + H{u) — H{u) 

Vt — Vt = DAm{v — ■!)) + F{u, v) — F{u, v) 

= G{u, v) - G{u, v) 

u — u = 0 
V — V = 0 


X € n, 

t > 0 

X G M, 

t > 0 

X G M, 

t > 0 

X G fl. 

t = 0 

X G M, 

t = 0 


Taking the dot product of the Vt — Vt equation with (y — v), and the Ut — itt equation with 
{u — u), and integrating over M and D respectively, yields 

- v\\Im + l|w - “lli.n) + £'||V(m - u)||i_n 

< ||f - v\\2M\\Fiu,v) - F{u,v)\\2,m + ll-a - u\\2,q\\H{u) - H{u)\\2y 
+ I|u - u\\2,m\\G{u,v) - G{u,v)\\2,m 

< K\\v - v\\ 2 ,M (||u - u\\ 2 ,M + Ik - ^IIz.m) 

+ K\\u - u\\ 2 .m (Ik - u\\ 2 ,m + Ik - v\\2,m) + +K\\u - 
<K{\\v-v\\jM + \\u-u\\jM) 

+ 2K\\u - u||2.M|k - v\\ 2 ,m + K\\u - 'u|| 2 ,n 

< 2K{\\v — v\\l^M + Ik ~ ^II 2 ,m) + Fi\\u — u\\ly 

From Lemma 2.4, for p = 2 and e = ~ have 

(6-4) Ik - w|l2.M < ^^11 V(w - u)||2,n + Cclk - u\\l^Q 

Using (6.4) 

~ ~ ~ + K{1 + 2Ce)\\u - u\\lf^ 

< C',,fc (Ik - f)||2,M + Ik - wllla) 

Observe, {u — ii) = {v — v) = 0 at t = 0 and (|k — u|l 2 a + Ik ~ '^112 m) — 0- Therefore, apply¬ 
ing Gronwall’s inequality, v = v and u = u. Hence system (3.1) has the unique global solution. □ 


Proof of Theorem 3.2: Recall that uq € VFp (D) and vq € Wp(M) with p > n, and uq^vq 
satisfies the compatibility condition for p > 3. From Sobolev imbedding (see [11], [19]), 
are bounded functions. Therefore there exists f > 0 such that |ko(’)||oo.n < r, |ko(’)lloo,M < f. 

For each r > f, we define cut off functions (jjr S (^“(IR^', [0,1]) and ijjr € (^(["((R^' x 
R’"),[0,1]) such that 4)r{z) = 1 for all k| < r, and 4’r{z) = 0 for all k| > 2r. Similarly 
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tpr{z,w) = 1 when \z\ < r and |w| < r, and ipr{z,w) = 0 when \z\ > 2r, or |w| > 2r. 
In addition, we define Hr = H(j)r,Fr = Fipr and Gr = Gipr- From construction, Hr{z) = 
H{z),Fr{z,w) = F{z,w) and Gr{z,w) = G{z,w) when \z\ < r and |w;| < r. Also, there exists 
Mr > 0 such that Hr, Gr and Fr are Lipschitz functions with Lipschitz coefficient Mr- Consider 
the “restricted” system 


(6.5) 


Ut = DAu + Hr{u) 

Vt = DAmV + Fr{u, v) 

= Gr{u,v) 

U = Uq 
V = Vo 


X G ft, 

t > 0 

X G M, 

t > 0 

X G M, 

t > 0 

X G ft, 

t = 0 

X G M, 

t = 0 


From Theorem 6.1, (6.5) has a unique global solution (ur, Ur). If ||u(-, t)||oo,n, ||i^(‘j i)||oo,M < 
r for all t > 0, then {ur,Vr) is a global solution to (3.1). If not, there exists Tr > 0 such that 

||Mr'(-,<)||oo,a + ||'yr(-,t)||oo.M <r Vt G [0, T^] 

and for all t > Tr there exists t such that Tr < t < t, and x G fl and z G M, such that 


\Ur{x,t)\ + \Vr{z,t)\ > T 

Note that Tr is increasing with respect to r. Let Tinax = limr->.oo Tr- Now we define {u,v) as 
follows. Given 0 < t < Tmax, there exists r > 0 such that t < Tr < Tmax. For all a; G ft, 
u{x,t) = Ur{x,t), and for all x G M, v(x,t) = Vr{x,t). Furthermore {u,v) solves (3.1) with 
T = Tmax- Also, uniqueness of (ur,Vr) implies uniqueness of (it,u). It remains to show that 
the solution of (3.1) is maximal and if T^ax < oo then 

limsup ||u(-,t)||oo,n + limsup ||u(-, t)||oo.M = oo. 

Suppose Tmax < cx) and set, 


limsup |lu(-,t)||oo.n + limsup ||u(-, t)||oo,M = R- 

^ ^^max ^ ^^max 


If i? = OO then (u, v) is a maximal solution. If i? < oo there exists L > 0 such that 

ll'«lloo.nx(O.T„ax) + l|T||oo,Mx(0.r„ax) - 
As a result, T 2 L > Tmax, contradicting the construction of T 2 l- Q 


Now we prove that under some extra assumptions that the solution to (3.1) is component¬ 
wise nonnegative. Consider the system 


Ut = DAu + H{u'^) 

Vt = DAmv + F{u'^,v~^) 
Fin 

(6.6) £>—= G(u+,i;+) 

u = Uo 


X gVL, 0 <t <T 
X G M, 0 <t <T 

X G M, 0 <t <T 

X gVl, t = 0 
X G M, t = 0 


V = Vo 
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where it’*' = max(u, 0) and u = — min(u, 0). 


Proposition 6.2. Suppose F^G and H are locally Lipschitz, quasi positive functions, and 
uq,vq are componentwise nonnegative functions. Then (6.6) has a unique componentwise non¬ 
negative solution. 


Proof. Note that P'(m+,p+), G{u^,v~^) and H{u~^) are locally Lipschitz functions of u and 
V. Therefore from Theorem 3.2 there exists a unique maximal solution to (6.6) on (OjTmax)- 
Consider (6.6) componentwise. Multiply the Vif equation by v~ and the equation by uj, 

dv' " 

(6.7) v~ = diV~ AMVi + v~ Fi(u+,w+) 

du ■ 

(6.8) uJ = djuj Auj + uJ Hj (u+) 

Since = 

^ ~ F,(u+,u+) 

- djUj Auj - uJHj{u+) 

Integrating (6.7) and (6.8) over M and Q, respectively, gives 


1 d 


1 d 


2^11^* (•>0ll2.M + ij:^l|wj + / jVu, I dcr + dj l \Vuj\ dx 


2dt" ^ 


f M 


= - [ u. dx - [ u, Gj{u~'~da - [ ,v~'~) 

Jn Jm Jm 


da 


Since F, G and H are quasi-positive and di,dj > 0, 






Therefore, the solution (u, v) is componentwise nonnegative. □ 


Corollary 6.3. Suppose F,G and F[ are locally Lipschitz, quasi positive functions, and 
uo,vo are componentwise nonnegative functions. Then the unique solution (u,v) of (3.1) is 
componentwise nonnegative. 


Proof. From Theorem 3.2 and Proposition 6.2, there exists a unique, componentwise non¬ 
negative and maximal solution {u, v) to (6.6). In fact {u, v) also solves (3.1). The result follows. 

□ 


6.2. Bootstrapping Strategy. The following system will play a central role in duality 
arguments. 


- - d 

(6.9a) ^' = 0 


{x, t) e M X (t, t) 
X G M, t = T 


{pt = —dAip — d 

.dip 

Kid^r- + K2T = W 
or] 


ip = 0 


{x, t) G it X (r, T) 
{x,t) G M X (t, T) 
X G fl, t = T 


(6.9b) 
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Here, p>l, 0<T<T,-d€ Lp{M x (r, T)) and i? > 0, and d G Lp(H x (r, T)) and -d > 0. 
Also d > 0, d > 0, and ki, /C 2 G M such that ki > 0 and KiK2 ^ 0. Lemmas 6.4 to 6.8 provide 
helpful estimates. 

Lemma 6.4. (6.9a) has a unique nonnegative solution 4' G x (r, T)) and there 

exists Cp^T > 0 independent of d such that 


Proof. The result follows from Theorem 3.5 and the comparison principle. □ 

Lemma 6.5. Let p > 1, ni > 0 and if ki = 0 then K 2 > 0. Suppose 4^ is the unique nonneg¬ 
ative solution of (6.9a). Then (6.96) has a unique nonnegative solution ip G Wp'^{n x (r,T)). 
Moreover, there exists Cp^T > 0 independent of d and d and dependent on d, d, ki and K 2 sueh 
that 


ll'dllp^(nx(r.T)) - C'p,T(||d||p_Mx(r,T) + II X (t.T) ) 


Proof The result follows from Lemma 6.4, Sobolev embedding and similar arguments of 
proof on page 342, section 9 of chapter 4 in [19], and the comparison principle. □ 


Remark 6. If p > n -\-2 and ki > 0, then is Holder continuous in x and t. See the 
Corollary after Theorem 9.1, (page 342) chapter 4 of [19]. 


Lemma 6.6. Suppose I > 0 is a non integral number, > 0, d > 0, d G (12 x [t,T]), 
d G C'’5(M X [t,T]), (p{x,T) G C^+\n) and 4- G C"+b^(M x [t,T]). Then (6.96) has a 
unique solution in x [r, T]). Moreover there exists c > 0 independent ofand d 

such that 


w\ 


{ 1 + 2 ) 

ax[r,T] 


< c 



(0 

Oxp.T] 


+ |4/| 


d+i) ^ 

Mx{t,T) J 


Proof. See Theorem 5.3 in chapter 4 of [19]. □ 

Lemma 6.7. Suppose 1 < p < oo, ki > 0, and r,s are positive integers. If q > p and 
2 — 2r — s — (n + 2) > 0 then there exists K > 0 depending on 12, r, s, n,p such that 

||7d[id«^||,,nx(.,T)<i^||7>llSx(..T) 

for all p G Wp'^{Cl X (r,T)). 

Proof. See Lemma 3.3 in chapter 2 of [19]. □ 


2r 


Lemma 6.8. 
f s < 2m - 

DID: 


Suppose 1 < p < oo, Ki > 0, and r,s,m are positive integers satisfying 
There exists c > 0 independent of p £ Wp^^’™‘{fl x (r, T)) such that 


G W, 


2m—2r—s — 


p (12) and \\p || 


(2m-2r--s-|) 

p,Q 


< c\\p lll,^. 




p,f2x(T,T) 


28 


V. SHARMA AND J. MORGAN 


In addition, when 2r + s < 2m — 


2m—2r — s— 


, ,, ,,(2m—2r — s—-) 

ll‘^ llp.Mx(r.T) " 

Proof. See Lemma 3.4 in chapter 2 of [19]. □ 


h X (r,r)) 

llpS(r.T) 


Lemma 6.9. Let p > 1, ki = 0 and suppose 0 < i? G Lp{n x (t,T)), and 41 is a unique 

solution of (6.9a). Then G Wp (M x {t,T)), and (6.9b) has a unique solution (p G 

Wp’^(fl X (r, T)). Moreover, there exists Opj- > 0 independent of tI and dependent on d, and 
K 2 such that 


ll‘^llp,ax(r.T) — ^P.r(ll'*^llp.nx(r,T) + ll^llp.M X (t,T) ) 


Proof. The result follows from Theorem 9.1 in chapter 4 of [19], Lemma 6.4, and Sobolev 
embedding. □ 

Remark 7. Ifp > , ki = 0 and ip satisfies system (6.9b), then p is a Holder continuous 

function in x and t. See the Corollary after Theorem 9.1, chapter 4 of [19]. 

Remark 8. By Lemma 6.4, Lemma 6.5, Lemma 6.8, and Sobolev embedding, we have 
2_^ 2_ 2. 

p(-,t) G Wp ’’ (LI), 4'(-,t) G Wp ’’ (M), and there exists c > 0 independent of p, 4/ such that 

II1’(g'’’) llp.n ^ — '^(ll'*^llp.nx(r.T) + ll^llp.Mx(r,T)) 

W^i-W) lliV^<c||^IUx(..T) 

respectively. Moreover, if p > n there exists c > 0 independent of p, 4^ such that 

ll<i’lloo.nx(r.T) < 


respectively. 


oo,M X (t,T) 


<c||4/(.,r)l| 


(2-i) 

p,M 


Lemma 6.10. Let 1 < p < n + 2 and 1 < q < There exists a constant C > 0 

depending on p,T — t, M and n such that if p € Wp "^(Ll x (r, T)), then 


dp 

dp 


q,Mx(r,T) 


<^11^11 


( 2 ) 

p,Qx{r,T) 


Proof. It suffices to consider the case when p is smooth in x [r, T], as such functions are 
dense in Wp’^(Q x (t,T)). M is a n — 1 dimensional manifold (p > 0). Therefore, for 

every f G M there exists > 0, an open set V C M" containing 0, and a diffeomorphism 

f) :V ^ B(l, Cj) such that '0(0) = 'i[({x G V : Xn > 0}) = B(^, e^) fl LI and 'i[’({x gV : Xn = 

0}) = R(^,e|) n M. Since 0 is a diffeomorphism, (tj}~^)n, the nth component of 0“^, is 
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differentiable in B{^, e^), and by definition of '0, (0 ^)ra(C) = 0 if and only if ^ S e^) fl M. 

Further, V(0“^)„(^) is nonzero and orthogonal to B{^,e^) fl M at each ^ G i?(^, e^) fl M. 
Without loss of generality, we assume the outward unit normal is given by 




V(r^)n(C) 

l(V0-i)4e)l 




We know, 


t) = t) ■ r;(0 V (e, t) G B{i e^) H PI x (r, T). 

Now in order to transform back to M", pick L > 0, such that 

E = [—L,L] X [—L,L] X ... X [—L,L] x[0,L] C V, and define (p such that 

'-V-' 

(n — 1) times 


ip{x,t) = - y^ip{ilj{x',z),t)'^D{ijj{x',z))r]{'il:{x',z)) dz W x = {x', z) £ E 

Jo 

where x' G [—L,L] x [—L,L] x ... x [—L,L]. We know (p G x (r,T)). Therefore from 

(n — 1) times 

Lemma 6.7, there exists 0 < a < L and > 0, depending on ff, n,p such that 


(6.9) 


f 

d(p{{x',a),t) 

ISc 

dXn 


dadt < ifj||¥^||J,'L(x.T) V 1 < r < 


(n + 2)p 


where Sa = E\x„=a x {,‘T,T) and Sx„ = i?|o<x„<a x (r, T). Using the fundamental theorem of 
calculus. 


1 

JEx{t,T) 

d<p{{x\0),t) 

dXn 

da dt < ( 
JSc 

d<p{{x',a),t) 

dXn 





d(p{{x', s 
dXn 


da dt 

d'^ip{{x',s),t) 


dxl 


da dt 


Using (6.9), 


f 

d(p{{x',0),t) ‘‘ 

'ex{t,T) 

dxji 


da dt < A^jdlv^llSxU.T)) 


[ 

d<p{{x',s),t) 

q-1 

d'^<p{{x',s),t) 


dXn 


dxl 


da dt 


Applying Holder inequality. 


[ 

dp{{x',0),t) 

Jex{t,T) 

dXn 


da dt < 


d(p{{x\s),t) 


dXr, 


(g-i)p 

p-i 


d(j dt 


d'^ipiix\s),t) 


dxl 


P 

da dt 
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Recall 



€ Lp(Sx„). So using Lemma 6.7 we have 


( 6 . 10 ) 


f 

dp{{x',0),t) 

Iex{t,T) 

dxji 


do- dt < ^(||</S||p^L(r,T))'^ 


Now, M is a compact manifold. Therefore there exists set A = {Pi,Pj^} C M such 
that M C Ui<i<NB{Pi,ep^). Let Vi, Ki and ai be the open sets and constants respectively 
obtained above when ^ = Pi. Then, 



PiGA 


Therefore, for some C > 0, depending only upon p, r, T, M and n, we get 


dp 

dr] 


q,Mx(T,T) 


< CM 


( 2 ) 

p,nx (r,T) 


for all 1 < g < 


(n + l)p 
n + 2 — p 


□ 

The following Lemma plays a key role in bootstrapping Lp estimates of solutions to (3.1). 


Lemma 6.11. Assume the hypothesis of Corollary 6.3, and suppose {u,v) is the unique, 
maximal nonnegative solution to (3.1) and T^ax < oo- If ^ j k and 1 < i < m, such that 
(RijT) holds, then there exists R'rmax > 0 such that 

\\uj{-,t)\\i^n + \\vi{-,t)\\i^M + l|Mj||i,Mx(0.T,„ax) < -f^T„ax for all 0<t< T^ax- 


Proof. For simplicity, take cr = 


( 6 . 11 ) 


pt = -dAp 
,dp 

d— = ap + l 
dr] 


1 in (14jl)- Let 0 < T < T^ax, and consider the system 
{x,t) G 11 X (o,r) 

{x, t) G M X (0, T) 


P = Pt 


X G fl, t = T 


where a is given in {Vijl), d > 0, and px G C'^+^(0) for some T > 0, is nonnegative and 
satisfies the compatibility condition 



apx + 1 on M X {T} 


From Lemma 6.6, p G C^+^’i +2 (Jl x [0,T]) and therefore by standard sequential argument 
p G X [0,r]). Also, note that g(s) = as + 1 satisfies g(0) > 0. Therefore, 

Proposition 6.2 implies > 0. Now having enough regularity for p on M x [0,T], consider 


Amp 


1 

■^det g 


d] {g''^ \Jdet g dip) 
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where g is the metric on M and g'^’^ is ith row and jth column entry of the inverse of matrix 
associated to metric g. Further let d = —ipt — dlS.M'-P- Then, 


Vi'd = 


/o JM 


'0 JQ, 

pT r 


/o Jo. 


Uj{—(pt — dA(p) + 
(p{uj^ - dAuj) + 


Vi{-ipt - HAm^p) 


'0 Jm 

'T r 


/o Jm 


ip{vit - dAMVi) - d 


'0 J M 


, 

Ut — + d 

dg 


/O JM 


9uj 

dg 




+ / 

Jo 

Using + 1 

/7 

/o Jm 


(x,0)v3(x,0) + f Vi{C,0)(f{x,0) - ( Uj{x,T)ipT - [ Vi{C,T)ipT 
Jm Jo Jm 


Uj < 


(pHj (u) ■ 


10 Jo 


10 Jm 


{F^{u,v) + Gj{u,v))ip 


Uj{x,Q)(p{x,0) + Vi{C,0)(p{x,0)- 


IM 


IQ JM 


Using (Ujl), 

( 6 . 12 ) [ [ Uj < f [ I3ip{uj + 1 ) + / f a{vi + l)ip 

Jo Jm Jo Jo Jo Jm 


lo 


Uj{x,0)(p{x,0) + Vi(C,0)ip(x,0))- 


IM 


lo Jm 


Now, integrating the Uj equation over U and the Vi equation over M, 
d 


dt 


I' 


Uj+ Vi]=d Auj + / Hj{u)+d I Avi+ / Fj{u,v) 


IM 


IM 


(6.13) 


< 13 [ (mj + 1)+ [ {Gj{u,v) + Fi{u,v)) 
Jo Jm 

<13 / {uj + 1) + a (uj + rij + 1) 

Jo Jm 


Integrating (6.13) over (0,t) with 0 <t <T < Tmax, and using (6.12), gives 

(6.14) f Uj{x,t) + f Vi{C,t) < P f f Uj + d f ( Vi+L{t) 

Jo Jm Jo Jo Jo Jm 

where 

L{t) = a\M\t + I3\n\t + a/3||v5||i.nx(o.t) + a^ll'^lli.Mx(o.t) + a\\uj{x, 0)||i,n • \\ip{x, 0)||oo,n 
+ lki(C,0)||i,M + a||Ti(C,0)||i,M • ||‘/?(a;,0)||oo.M + ||uj(x, 0)||i,n 

d{t) = a^||(/?||oo.Mx(o.t) +a + a!||0||oo.Mx(o.t) and /3(t) = 13 + a^||(/?||oo,ax(o.t) 

Applying Generalized Gronwall’s inequality to (6.14) gives the bound for the first two integrals 
on the RHS of (6.14), and then substituting this bound gives 

f Uj{x,t) + [ Vi{C,t) < L(t) + f (a(s) +/3(s))L(s) exp f d{r) + P{r)dr\ ds 
Jo Jm Jo \Je / 

< Gt 

— J- max 
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for all 0 < i < T < T^ax- Substituting this estimate of Uj on fi and Vi on M in (6.12) yields 


n Uj < ^(||</3|| C50, nx(o.T)l|wjl|i,ax(o,T) + l^l'r||(/3||oo,ax(o,T)) 

+ O' (||‘/5||oo.Mx(0,T)|ki||l,Mx(0,T) + |■^^|T"||</^||oo,Mx(0,T)) 

+ l|Mj(-,0)||i,a||(/5(-,0)||oo,n + ||ui(-,0)||i,M||<p(-,0)||oo.M + ||ui||i,M||0||oo,M 
Since T < Tmax is arbitrary, the conclusion of the theorem holds. □ 


Lemma 6.12. Assume the hypothesis of Corollary 6.3 holds. Suppose (u,u) is the unique, 
maximal nonnegative solution to (3.1) and T^ax < oo. If ^ < j < k and \ < i < m, such that 
and {Vij2) holds, and for q > 1, Vi £ Lq{M x (0,Tinax)), then Uj G Lq{M x (0,Tinax)) 
and Uj £ Ijq{ft X (0,X'rnax)). 

Proof. Let 0 < t < T < Tmax. Multiplying the Uj^ equation by we get 


q—1 

^ = 


= d 


Jo Jo. Jo Jo. 

_i duj 


0 Jm 




0 JQ. 


u] ^Hj(u) 


Using (Vij2) 


f^+df / < iLg / / u] ^{uj+Vi + l)+P f f iuj+l)u'j ^ 

Jo, <7 Jo JQ. H Jo Jm Jo Jn 


/ )Uo 

In 9 

ft 


<Kq 


0 JM 


u‘j + ViUj ^ +Uj M 


Q I 9—1 

U j + M j 


(6.15) 

Applying Young’s inequality in (6.15) 


In 9 


/ ^ + d 
In 9 


4(9-1),^ f 


0 Jn q 


\Vu^A^<Kq 


iq-2 


u® + f/3 + t|U| — 


0 JM 


(6.16) 


tn 9 


^+Kq[- 


qj Jo Jn 

,t\M\ 


qj Jo Jm 


Also, for 1 < g < oo, for all e > 0 and t < T < T^ax, from Lemma 2.4, for v = there exists 
Ce > 0 such that. 


(6.17) 


< C, / u’^q+e 
lo JM Jo Jn 


r‘ 

to Jn 


|Vu| 


Using (6.17) and (6.16) for appropriate e > 0, gives 


(6.18) 


l_d '■* 
q dt 


0 Jn 


u]<K, 


0 Jn 


u] + K2{T) 
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for 


K2{T)=Kg 


T\M\ 


.Q/ Jo Jm Q 

and Ki > 0 depending on t, where t < T < T^ax- Therefore from Gronwall’s Inequality 


(6.19) 


Uj'^{x,t) < K 2 {T) + / Ad(s)A" 2 (s) exp / ki{r)dr\da 


To obtain estimates on boundary, we use (6.16) to obtain 

pT p / 2 ~ 

r n ( q 


'0 Jo. 




( 6 . 20 ) 




+ 


ul + e 


0 JM 


^ \4.d{q — 1) 


P + T\n\ 


Ad{q-l)J \ Jn q 


' eKg{- 


Using (6.17), (6.20) and (6.19) we have. 


/o Jm 


u] < C. 


10 Jo 


u] + ■iKg 


qj Jo Jm 


f q^ 


Ad{q-l)J Jo Jm ^ \Ad{q-\) 

T\M\\ 


qj Jo Jo 
T\M\\ 


p + T\nj- 


qj Jo Jo 


AM 1 W —+^Kg\- 

Ad{q-l)J \ Ja q 


Now choosing e such that 


l-3iC„ 


Ad{q - 1) 


qJ Jo Jm 


e > 0 


and using the estimate above for Uj on (11 x (0,T)), we have Uj G Lq{M x (0,T)). Since T is 
arbitrary, Uj G Lq{M x {Q.Tmax)) □ 


Lemma 6.13. Assume the hypothesis of Corollary 6.3, and suppose (u, u) is the unique, 
maximal nonnegative solution to (3.1) and < oo. If A < j < k and 1 < i < m so that 

(Vijl) and (Vij2) hold, then for all p > 1 and 0 < T < Tnia,x, there exists Cp^T > 0, such that 

l|Wj||p,nx(0,T,„ax) + ll^*llp,Mx(0,T,„ax) 

< Cp,Tmcx (ll'“jlll.Mx(0,T,„ax) + II l.nx(0,T,„ax) + 11 111 .M X (O.T,„„x ) ) 


Proof. First we show there exists r > 1 such that if g > 1 such that Uj G Lq{rt x (0, Tmax)) 
and Vi G Lqi^AI X {0,Pmax)) then Uj G Lpq{yi X (0,Pmax}} and Vi G Lj-qi^lM X (0,Pmax}}- 
Consider the system (6.9a) and (6.9b) with = 0, K 2 = 1, > 0, d G Lp(M x (0,Tmax)) 

with ll")? llpyMx(0,Tmax)) ~ ^ S Lp(fl X (9,Tmax)) with ||l? lip,(f2x(0,Tmax)) ~ 

Multiplying Uj with d and Vi with d and for 0 < T < Tmax, integrating over 11 x (0,T) and 












34 


V. SHARMA AND J. MORGAN 


M X (0, T) respectively, gives 





Uj{-ipt-dA(p)+ / 


lo Jn 

rT r 


lo Jm 

rT r 


lo Jn 


tp{uj^-dAuj) + / - dAMV^) 


lo JM 


Jo Jm or] Jq jfj 07 ] Jq 

+ f Vi{x,0)'^{x,0) - [ Vi(x,T)'^{x,T) - f Uj{x,T)(p{x,T) 
Jm Jm Jn 


Since '^{x,T) = 0 and (p{x,T) = 0, 


/7 


Ujd- 


fl 

Jo J M 


k'O < 


lo Jo 


-d 


ipHj {u) + 


{Fj{u,v) + Gi{u,v))'i> 


lo JM 


Jo JM or] Jq 

+ [ Vi{x,0)'il{x,0) 
JM 


Using {Vijl), 


( 6 . 21 ) 




I3ip{uj + 1) + 



a{uj +Vi + 1)'!' 


-df [ Uj^ + [ Uj{x, 0 )ip{x,Q) 

Jo Jm 07] Jq 

+ [ r;i(a;,0)'i;'(a;, 0) 

JM 


Now we break the argument in two cases. 

Case 1: Suppose q = 1. Then Uj G Ti(U x (0,Tj„aa;)) and Uj,Vi € Li{M x (0,Tmaa:))- Let 
e > 0 and set p = n + 2 + e. Set p' = (conjugate of p). Remarks 6 and 8, and Lemma 

6.11 imply all of the integrals on the right hand side of (6.21) are finite. Application of Holder’s 
inequality in (6.21), yields Vi € Lpi{M x (0,T)), and there exists Cp^T > 0 such that 


l|Millp',ax(0.T) + ||w*||p'.Mx(0,r) < Cp^T{\\Uj\\i^nx(0,T^aA + lki||l,Mx(0,T„„,.) + Iki || l.Mx (0,T^a„)) 

Since T < T^ax is arbitrary, therefore, Lemma 6.12 implies Uj G Lp'{M x (0,Tmax))- So for 
this case, r = 

Case 2: Suppose g > 1 such that uj G Lg(H x (0,Tmaa;)) and Uj,Vi G Lq{M x (0,rmaa;))- 
Recallp > 1, 0 < d G Lp{M x (0, T^ax)) with ||d \\p,]Mx{o.Tma.^)) = 1 and 0 < -d G Lp{Q x (0, Tmax)) 
with ||d ||p,(nx(o,T„„,,)) = 1- P' = Note T < Tmax is arbitrary. Applying 

Holder’s inequality in (6.21) and using Lemma 6.10, yields 


Ullp'.nx(0.T™ 


) + ll^*llp',Mx(0,T„„x) 

< Cp,T,nax{\\ll3\\q,Ox{0,T^a.x) + I l^i 11 <?. M X (O ) + \\Uj\\q,Mx{,0,Tm^^)) 


provided p' < 777 ■ 


So, in this case, r = 


(n+2) 
n+1 * 
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Now, by repeating the above argument for rq instead of q, we get Vi € Lrmq{M x (0, Tmax)), 
Un € Lr^„(n X (0,Tjnax)), for all m > 1. As r > 1, lim r'^q —>• oo, and as a result, Vi G 

Lp{M X {Q,Tjnax)) for all p > 1. Hence from Lemma 6.12, Uj G Lp{M x {Q,Tmax)) and 
Uj G Lp(fi X (0,Tmax)) for all p > 1, and there exists Cp^r > 0 such that 

l|Mjllp.nx(0,T) + ll'yillp.Mx(O.T) < C'p.T (l|ujllg.Mx(O.T,„„x) + \Wj\\q,nx{0,T,r,ax) + lkillg.Mx(O.T™„x)) 

Again as T < Tmax is arbitrary, we get 

l|uillp,nx(o.T^„,,) + lki|lp,Mx(o.r^„,,) 

< Cp,T,ncx (ll'ailll,Mx(0,Tmax) + l|■allll,Ox(0,Tmax) + 11 111 ,M X (0,T„ax ) ) 


□ 

6.3. Global Existence. 


Proof of Theorem 3.3: From Theorem 3.2 and Corollary 6.3, we already have a compo¬ 
nentwise nonnegative, unique, maximal solution of (3.1). If T^ax = oo, then we are done. So, 
by way of contradiction assume Tmax < oo. From Lemma 6.13, we have Lp estimates for our 
solution for all p > 1, on H x (OjTmax) and M x (0,rmax)- We know from (14j2) and (14^3) 
that Fj and Gi are polynomially bounded above for each i and j. Let U and V solve 


( 6 . 22 ) 


Ut = djAU + ^{uj + 1) 

V) — diA][^V + Kj(^Uj Vi 1)^ 

dU 

U = Uo 
V = Vo 


{X,t) G O X [Q.Tmax) 
{X,t) G M X {0,Tmax) 

{x,t) G M X {0,Tmax) 

X G fl, t = 0 
X G M, t = 0 


Here, dj and di are the jth and ith column entry of diagonal matrix D and D respectively. 
Also, C/q and Vq satisfy the compatibility condition, are component-wise nonnegative functions, 
and {uo)j < Uq and (r!o)i < Vq. For all q > 1, Kf{uj + Vi + 1)* and Kg{uj -|- -|- 1) lie in 

Lq{M X {0,Tmax))- Using Theorem 3.6, the solution of (6.22) is sup norm bounded. Therefore, 
by the Maximum Principle [26], the solution of (3.1) is bounded for hnite time. Therefore 
Theorem 3.2 implies Tmax = oo. □ 

7. Examples and an Open Question. In this section we give some examples to support 
our theory. 


Example 1. As described in [35], during bacterial cytokinesis, a proteinaceous contractile, 
called the Z ring assembles in the cell middle. The Z ring moves to the membrane and contracts, 
when triggered, to form two identical daughter cells. Positiong the Z ring in the middle of the 
cell involves two independent processes, referred to as Min system inhibition and nucleoid 
occlusion ([32], [33] Sun and Margolin 2001). In this example, we only discuss the Min system 
inhibits process. The Min system involves proteins MinC, MinD and MinE ([28] Raskin and de 
Boer 1999). MinC inhibits Z ring assembly while the action of MinD and MinE serve to exclude 
MinC from the middle of cell region. This promotes the assembly of the Z ring at the middle of 
the cell. In [35] the authors considered the Min subsystem involving 6 chemical reactions and 
5 components, under specific rates and parameters and performed a numerical investigation 
using a finite volume method on a one dimensional mathematical model. Table 7.1 shows the 
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assumed chemical reactions. The model was developed in [35] within the context of a cylindrical 
cell consisting of 2 subsystems; one involving Min oscillations and the other involving FtsZ 
reactions. The Min subsystem consists of ATP-bound cytosolic MinD, ADP-bound cytosolic 
MinD, membrane-bound MinD, cytosolic MinE, and membrane bound MinD:MinE complex. 
Those are denoted and E : respectively. This essentially 

constitutes the one dimensional version of the problem. These Min proteins react with certain 
reaction rates that are illustrated in Table 7.1. These reactions lead to hve component model 


Table 1: Reactions and Reaction Rates 


Chemicals 

Reactions 

Reaction Rates 

Min D 

nADP toATP 

^cyt ^ ^cyt 

Rexc = kllD^yr] 

Min D 

ToATP nATP 

^cyt ' ^mem 

Rncyt = k2[D^Jin 


pATP pATP 

^cyt ' ^mem 

Romem = k3[D^l^][D^Jt^] 

Min E 

rp 1 t^ATP 77' . r^ATP 

^cyt + ^mem ^ ^ ■ ^mem 

Rpcyt = k4[Ecyt][D^^] 


, 7-iATP ^ ^ ^ t-iATP 

P^cyt + ^mem ^ ^ ■ E'mcm 

Rpmem = fcs [E,,*] [A : 

Min E 

IT' . T~\ATP 77' 1 T~)ADP 

h/ . -LfjYiQYn ' ^ ^cyt 

Rexp = ke[E : 


with {u,v) = {ui,U 2 ,U 3 ,vi,V 2 ), where 




D 


ATP' 

cvt 

ADPl 


U = \U2\ = \ [D^y^ 
\ [Ecyt 


,v = 


V 2 J [[b-.d: 


ATP] 
mem J 


D = ^ D 

V U ^E\Dmem, 


' ^Dcyt 0 0 

0 ^ADyct 0 

0 0 ^Ecyt y 


G{u,v) = I G 2 {u,v) 
.G3(u,v), 


Epcyt Rd mem 
Rexp 

\Rexp ^Ecyt Re mem y 


jpf \ f _ f Rpcyt Rd mem REcyt 

yF2{u,v)J ~ -Re.p + REcyt +R, 


—k2Ui — k3ViUi 

keV2 

\keV2 — k4U3Vi — k3ViU3V2^ 


— Re mem 
Re mem 


_ fk2Ui + k3ViUi — k4U3Vi — k5ViU3V2" 
\ —keV2 -f k4U3Vi -I- k3ViU3V2^ 


/Hi{u)\ / Rcxc \ ( kiU2 \ 

H{u) = H2(u) = -Rcxc = -kiU 2 , 

\H3{u)I \ 0 j \ 0 j 

and uq = {uoj) G Wp{fl), vq = (uoi) & Wp{M) are componentwise nonnegative functions with 
p > n. Also, uq and vg satisfy the compatibility condition 

duo 
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Here expressions of the form ka and ap are positive constants. Note F, G and H are quasi 
positive functions. In the multidimensional setting, the concentration densities satisfy the 
reaction-diffusion system given by 


Ut = 

= DAu -f 

H{u) 

X gVl, 

0 

< 

t 

< 

T 

Mt = 

= DAmv 

+ F{u,v) 

X G M, 

0 

< 

t 

< 

T 

Dp, 

Ol] 

= G{u, v) 


X G M, 

0 

< 

t 

< 

T 

u = 

= Ug 


X G fl. 

t 

= 

0 



V = 

= Mo 


X G M, 

t 

= 

0 




Our local existence result holds for any number of finite components. Therefore, from Theorem 
3.2, this system has a unique maximal componentwise nonnegative solution. In this example, 
if we take two specific components at a time, we are able to obtain Lp estimates for each of 
the components. For that purpose we apply our results to (m 3,U2), U 2 and (ui,ui). In order to 
prove global existence, we assume T^ax < oo. Otherwise, we are done. 

Consider (ug, U2). It is easy to see that for j = i and i = 2, the hypothesis of Lemma 6.13 is 
satisfied, since G3 -I- F2 < 0, G3 is linearly bounded, and H 3 =0. As a result, M3 G Lp{^T,nax) 
and V 2 G for all p > 1. Using Theorem 3.6 and the comparison principle, U 2 is 

Holder continuous on for p > n + 1. Finally, consider (mi,mi). Clearly for j = 1 and 

1 = 1, the hypothesis of Lemma 6.13 is satisfied, since Gi + < 0, Gi is linearly bounded, 

and Hi is bounded. Therefore, mi G and vi G Lp{MT„^^) for all p > 1. 

We already have for all 1 < i < 2 and 1 < j < 3, {uj,Vi) G Lp{^Tmax) ^ ^p{^^T,nax) 
for all p > 1. Therefore there exists p > 1 such that Gj G Lp{i2T^^^) for all p > p, and 
Fi G Lp[MT^^^) for all p> p. Consequently from Theorem 3.6, the solution is bounded, which 
contradicts the conclusion of Theorem 3.2. As a result, the system has a global solution. 
Example 2. Consider the model considered by Ratz and R6ger[29] for signaling networks. 
They formulated a mathematical model that couples reaction-diffusion in the inner volume to 
a reaction-diffusion system on the membrane via a flux condition. More specifically, consider 
the system (3.1) with k = l,m = 2, where 


G(m, v) 


h 1^1 I 


- Ml - M2)+ -f b_eV2, 


H{u) = 0 


F{u, v) 


fFi{u,v)\ 

\F2{u,v)) 



1 - 
1 - 


Ksvigo A 






Ksvigo 


and Mo = (uoj) G Wp‘^\n), vq = (moi) G W^\m) with p > n are componentwise nonnegative. 
Also, Mo and vg satisfy the compatibility condition 

D-^ = G{uo,vo) on M. 

Here ka, Ka, go,Cmax,b-G are same positive constants as described in [29]. We note F,G and 
H are quasi positive functions. From Theorem 3.2, this system has a unique componentwise 
nonnegative maximal solution. In order to get global existence, we assume T^nax < 00 . In 
order to obtain Lp estimates for each of the components, consider (m, M 2). It is easy to see that 
G + F 2 < kg, H = 0, and G is linearly bounded above. So the hypothesis of Lemma 6.13 is 
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satisfied. As a result, u G ^p(^T,r,ax) V2 G for all p > 1. Now vi, satisfies the 

hypothesis of Theorem 3.5. Therefore vi G fo^ p > 1. 

We already have for all 1 < f < 2, (u, Vi) G Lp{GtT^^^ ) x Lp{MT^^^ ) for all p > 1. Therefore 
G G Lp{ilTma.x) fo'' 8-11 P ^ 1; s-nd Pi £ PpiP^Tmax) fo^ ^ 1- Cousequeiitly, from Theorem 
3.6, the solution is bounded, which contradicts the conclusion of Theorem 3.2. As a result the 
system has a global solution. 

Example 3. We look at a simple model to illustrate an interesting open question. Consider 
the system 


(7.1) 


ut = Am 
Vt = I^mv + 
du 
drj 

u = Uq 


2 2 
= —u V 


V = Vo 


X Gft, 0 <t <T 
X G M, 0 <t <T 

X G M, 0 <t <T 

X G ft, t = 0 
X G M, t = 0 


where Uq and vq are nonnegative and smooth, and satisfy the compatibility condition. Clearly 
H(u) = 0, G{u,v) = and F{u,v) = —u^v"^ satisfy the hypothesis of Theorem 3.3 with 
F + G < 0 and G{u,v) < 0. Therefore (7.1) has a unique global componentwise nonnegative 
global solution. However, suppose we make a small change, and consider the system 


(7.2) 


Ut = Am 

Vt = AmV — M^M^ 


du 

drj 


2 2 
= U V 


U = Uq 


V =Vq 


X G ft, 

0 

< 

t 

< 

T 

xGM, 

0 

< 

t 

< 

T 

xGM, 

0 

< 

t 

< 

T 

X G ft, 

t 

= 

0 



X G M, 

t 

= 

0 




Then we can show there exists a unique maximal componentwise nonnegative solution. We 
can also obtain Li estimates for u and v. Furthermore, it is easy to see that v is uniformly 
bounded. But our theory cannot be used to determine whether (7.2) has a global solution, 
and this remains an open question. More generally, it is not known whether replacing Gj in 
condition (Vij2) with Ft will result in a theorem similar to Theorem 3.3. 
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